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Abstract 

A four-fermion model in 2+1 dimensions describing N Dirac fermions interacting via SU(A^) 
invariant N'^ — 1 four-fermion interactions is solved in the leading order of the expansion. 
The expansion corresponds to 't Hoofts topological 1/N expansion in which planar Feynman 
diagrams prevail. For the symmetric phase of this model, it is argued that the planar expansion 
corresponds to the ladder approximation. A truncated set of Schwinger-Dyson equations for the 
fermion propagator and composite boson propagator representing the relevant planar diagrams is 
solved analytically. The critical four-fermion coupling and various critical exponents are determined 
as functions of N. The universality class of this model turns out to be quite distinct from the Gross- 
Neveu model in the large limit. 
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I. INTRODUCTION 



Already in the early seventies K.G. Wilson H] argued that four-fermion theories in dimen- 
sions 2 < d < A are nontrivially renormalizable at least in leading order of the 1/N expansion, 
with the number of fermion flavors. The model studied by Wilson is nowadays referred 
to as the d dimensional generalization of Gross- Neveu model [Q] . In the Gross- Neveu model 
with fermion flavors there is one composite or auxiliary field a oc — J2a=i '4'a'4'a describing 
the scalar bound states; the number of fermions is taken to be much larger than the num- 
ber of light bound states. In the 1/A^ expansion for such models, Feynman diagrams with 
fermion loops dominate over other types of diagrams. 

The Gross-Neveu model in 2+1 dimensions is well-known to exhibit dynamically sym- 
metry breaking whenever the four-fermion coupling exceeds some critical value. About a 
decade ago renewed interest in the Gross-Neveu model led to a number of papers establishing 
the renormalizability of the Gross-Neveu model and related models in the 1/A^ expansion 
due to the presence of a finite ultraviolet stable fixed point and conformal invariance at 
this critical point 0, ||, |^, |, |^. Probably the most extensive study of the Gross-Neveu 
model, was performed by Hands, Kocic, and Kogut In their paper analytical results 
in next-to-leading order in 1/A^ were verified in numerical lattice simulations. The 1+1 
dimensional four-fermion models cannot exhibit dynamical symmetry breaking, due to the 
Coleman-Mermin- Wagner theorem, whereas 3 + 1 dimensional models usually suffer from 
the triviahty problem ||^, |T^. 

In this paper a SU(A^)x U(l) invariant four-fermion model, closely related to the Gross- 
Neveu model, is studied in 2+1 dimensions. In this model the four-fermionic potential 
contains A^^ — 1 terms, i.e., the four-fermion interactions are in the adjoint representation of 
the SU(A^) symmetry. Thus the model presented here describes the interaction of A^ Dirac 
fermions with A^^ — 1 scalar composite states. Consequently the large A^ treatment will be 
quite different from the Gross-Neveu model. Instead of the large A^ expansion a la Wilson, 't 
Hoofts topological 1/A^ expansion |Tl| will be adopted. We assume that the 1/A^ expansion 
of 't Hooft is particularly useful for theories with a global U(A^) symmetry containing fields 
with two U(A^) or SU(A^) indices, e.g., see Refs. 0, |T^, |TT], |T^. In such an expansion 
planar diagrams dominate over diagrams with topologies other than planar. 

The leading large A^ critical behavior of the proposed model is studied, with the main 
focus on the determination of the scaling behavior of the fermion wave function and of 
the A^^ — 1 scalar or a boson propagators. Once the anomalous dimensions of the fermion 
propagator and the a boson propagators are determined, the gap equation for the order 
parameter will be studied. The dynamical breaking of the SU(A^) symmetry to U(l) gives 
rise to the Goldstone realization with A^^ — 2 Nambu-Goldstone bosons and one massive 
scalar in the broken phase. 

The main motivation for studying this particular three dimensional four-fermionic model 
is it possible relationship with low temperature behavior of the 2+1 dimensional Hubbard- 



Heisenberg and t — J models near the antiferromagnetic wave vector [[13| , |I4[ [T^. 

The setup of the paper is the following. The model Lagrangian is introduced in the next 
section. We present a truncation scheme for the Schwinger-Dyson (SD) equations which gen- 
erates all planar and thus leading 1/A^ diagrams in Sec. |ITl|. This closed set of SD equations 
is subsequently solved analytically in Sec. fV[ In Sec. the dynamical symmetry breaking 
is analyzed and the critical exponents, hyperscaling, and universality are discussed. Next- 
to- leading order 1/A^ corrections are discussed in Sec. |VJ The conclusions are presented in 
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Sec. [VI]| . In the Appendix , the apphcabihty of the topological expansion is motivated 
for this 2+1 dimensional model. 



II. THE MODEL 



We will consider the following 2+1 dimensional four-fermion model given by the La- 
grangian 

G 

C = Ipjdlpa + 17 II (1) 
^ A=l 

where the index a = 1, . . . ,N labels the N fermion flavors. The fermions are described by 
four- component Dirac spinors. The Lagrangian Eq. (|^) is globally invariant under SU(iV) x 
U(l), see also Ref. 0. By the way, in Refs. p|, |10[ a U(A^)x U(A^) invariant version of 



Eq. (|^) is described in 3 + 1 dimensions. The Hermitian generators, r"^, of the SU(iV) 
symmetry satisfy 

Tr[r^] = 0, TT[r\^] = 5^^ = ^o.sS,, - -S^,S,s. (2) 

A=l 

After a Hubbard-Stratonovich transformation, the Lagrangian (|I|) reads 

C = i^JdiJo, - i'ar^pi^/3(TA - 7^ H (3) 
A A 

with the real auxiliary bosonic fields = —GipT^ip conveniently describing the composite 
fermion-antifermionic degrees of freedom. Another equivalent representation with aij = 
J2a '^^^ij can be adopted [l^ , see also the Appendix . Along the lines of Ref. [|l^ , A^^ — 



1 cr boson propagators A^(p) and Yukawa vertices T'^{p + q,p) ("fully amputated") are 
introduced. The bare Yukawa vertices are •y^ = t^I. Since in the symmetric phase the 
propagators are identical, the index A is often discarded. 

The dynamical symmetry breaking SU(A^) — ^ U(l) is described by the order parameter 
((T3) ~ —G{ipT^ip) , which is a particular choice for the symmetry breaking term. In principle 
we have the freedom to choose from N"^ — 1 order parameters ipT^ip. In order to "force" the 
symmetry breaking in the direction a bare mass term —hipr^ip is added to the Lagrangian 
in Eq. (|^). In the end the limit — >■ can be taken. 

In the fermion propagator the symmetry breaking is given by the mass function M which 
is introduced as follows 



A{p)P5o,p + VNMtIp 

The function A{p) is referred to as the fermion wave function. The fermion dynamical 
mass TTLdyn is defined as the mass pole in the fermion propagator; it is related to M via 
^^{f^dyn)fTi\n ~ Both the fcrmiou dynamical mass nidyn or the mass pole mo- in the a 
boson propagator can be considered as the appropriate inverse correlation length character- 
izing the infrared mass scale of this model. 



3 




FIG. 1: The planar approximation; the dashed hnes are the real a boson propagators, the black 
lines are fermion propagators. The blobs correspond to the fully dressed propagators. The dark 
blob in the fermion propagator corresponds to a wave function part being proportional to j5, whereas 
the white blob with a cross corresponds to a mass insertion proportional to the generator. 



III. LARGE N PLANAR AND BARE VERTEX APPROXIMATION 

The present model describes N species of fermions interacting via A^^ — 1 four-fermionic 
interactions in the so-called adjoint representation of the SU(A^) flavor symmetry. There- 
fore, when N is large, we assume the applicability of 't Hoofts topological expansion 
111 , see also Ref. [|1^. Moreover, since the U(iV) symmetry is global, we assume that 



the corresponding topological expansion is independent of the space-time dimension d. 
The applicability of 't Hoofts expansion to the present model is motivated in the Ap- 
pendix . The topological expansion states that the Feynman diagrams can be classified 
in terms of surfaces with a specific topology. The so-called planar Feynman diagrams, with 
the fermions at the boundary, describe the leading contributions to the Green functions. 
Diagrams with other than planar topological structures are multiplied by factors of at least 
and in the limit of large A^, their contribution can be neglected with respect to planar 
graphs. For instance diagrams with one handle on their representing surface are suppressed 
by a factor relatively to their planar counter parts, whereas diagrams with a hole are 
suppressed by a factor 

In this paper we are mainly interested in the SD equations for the a boson propagator 
Act and the fermion wave function A. In the Appendix , it is shown that the leading 
large behavior for Ao- is given by those diagrams which do not contain Yukawa vertex 
(Fct) corrections, see the first two truncated SD equations in Fig. |l]. In other words, only 
self energy corrections (S) and a boson vacuum polarization corrections (IIo-) should be 
considered. One can straightforwardly verify that the bare vertex approximation (see Fig. 
r^(p + ~ T"^!, generates only planar and thus leading large graphs for A^- and A. 
Here we mention that the propagators of the real scalar fields a"^ are represented by dashed 
lines, whereas the propagators of the Hermitian fields a^p of the Appendix are represented 
by double lines. 

For the symmetry breaking fermion mass M only the tad-pole graph contributes in the 
leading large A^ limit, see Fig. |T[ At a first glance this might seem inconsistent, since parts 
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FIG. 2: The bare vertex approximation. 

of the fermion propagator are now described by different truncations. However, tlie self- 
consistency of the truncated set of SDEs given in Fig. |l] can be verified using the SU(A^) 
Ward identities, e.g., see Ref. |]10|- These Ward identities relate the iV^ — 1 Yukawa vertices 



V^{p,p) at zero a boson momentum transfer to the fermion dynamical mass M{p) in the 
following way: 



r^(p,p) = r^^, A = l,...N\ A ^3, (5) 



rj,..) - (a) 

when the symmetry breaking is in the direction. Thus if we apply the bare vertex 
approximation = (Fig. we get the equation M = (0-3). This equation is indeed 
represented by the third SD equation in Fig. 0. 

Summarizing; the leading large N behavior can be described consistently in terms of the 
truncated set of SDEs given in Fig. |l|. 



IV. SCALING BEHAVIOR AT CRITICALITY 

Having introduced the leading large planar or bare vertex approximation for this model, 
we proceed by studying the truncated set given in Fig. |l| in the symmetric phase first. These 
equations are in Euclidean formulation 

^-\p) = -^ + n.(p), (7) 

r d?k 4(A;2 + k ■ p) , . 

"^^^ ~ Je (2^ {k + pYk^A{k + p)A{k) ' ^ ' 

and 

A{p) = 1 - NE4p), (9) 

For the dimensionless coupling constant g = 2GA/tt'^ there exists an ultraviolet stable fixed 
point gc above which the SU(A^) is dynamically broken. Exactly at the critical coupling 
g = gc the infrared mass scale vanishes and the above propagators will be given in terms of 
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pure power-laws. Therefore, we first consider the two coupled SD equations at the critical 
point {g = Qc), where we can safely assume a power-law form for the fermion wave function 



Aip) = (p7A 



2/A2^-C/2 



(11) 



where A is the ultraviolet cutoff, and where both the critical coupling Qc and the exponent ( 
are to be determined from the above SD equations and will depend on N via Eq. (|). Using 
the above ansatz for the fermion wave function, we subsequently solve the SD for the a 
boson propagator. For the model to be nontrivially renormalizable, the a boson propagator 
should also exhibit a pure power-law form at the critical coupling; 



^Ap) oc ( - 



2-ry 



(12) 



where t] is the anomalous dimension of the a boson propagator. Therefore, in what follows, 
we solve the SD equations ([^)-([T0|) assuming the power-law forms ( PD and (|12D in order to 
obtain the exponents ( and 77 and the critical coupling Qc as function of A^. 
Inserting Eq. (pUD in Eq. (|]) we obtain the integral 



n.(p) 



2A-2C 



dk 



dQ 



{k^ + k- p) 



471 A;-C(A;2+p2 + 2fc-p) 1-^/2' 



The angular integral is straightforward 
na(p) = — — I dk 



2C(2 + C)A;p 



+ Ck^ + Ckp-p^){k + pY 



- {k^ + (k^ - (kp - p^){\k - p\)^ 



where k = y/W and p = \fp^. We split this integral into three parts; 

n^(p) = (^) \j\ \m,x)^j2{c,,x)^j^{c,x)\ 



with X = pI k and where 
with 



(1 + 

2C(2 + C) 
1 

2(2 + 
1 

"2C(2 + C) 



i 1 

jdtt^^'h{t) + Jdtt-'<-'Mt) 



dtt^h{t)+ / dtt-^'^-^Ut) 



dtt^'^fiit) + / dtr^^-^fiit) 



LO 



(13) 



(14) 



(15) 



(16) 
(17) 
(18) 



h{t) = (1 + t)< - (1 - t)^ h{t) = (1 + tf + (1 - t)< 



(19) 
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The integrals can be performed straightforwardly using the Table of Integrals [|T7| and the 
J function can be expressed in terms of the F function and the hypergeometric function 
2-Fi = F . In this way, the function J\ can be expressed as (with x = p/ K) 



2C(2 + C) 



[F(-2 - 2C, -C; -1 - 2C; -x) - F{-2 - 2C, -C; -1 - 2C; x) 



(2 + 20 
(2TC)^ 

r(2 + c)r(i + c) r(-2-2C)r(i + c) 



+ (^F(2 + C, -C; 3 + C; -1) - (^^(-2 - 2C, -C; -i - 2C; -i) 



(20) 



r(3 + 2C) r(-i-c) ]■ 

Subsequently J^2{(,x) and J^siC^x) are computed in a similar manner, giving 

^2(C,^) = ^^^^|^^^[F(-i-2C,-C;-2C;-x) + F(-i-2C,-C;-2C;a;)] 

+ (TTO^^^ ^ ^' ^ + ^' "^^ " (rT2C)^^"^ " "^^ 
r(i + c)r(i + c) r(-i-2C)r(i + c) i . . 

r(2 + 2C) r(-c) J' ^ ^ 



and 



JsiCx) = -^^J-^|^[F(-2C,-C;l-2C;-x)-F(-2C,-C;l-2C;x) 

+ ^F{C, -C; 1 + C; -1) - ^^(-2C, -C; i - 2C; -i) 

_ r(C)r(i + c) r(-2C)r(i + c) l 
r(i + 2C) r(i-c) ]■ 



(22) 



The condition for the existence of these integrals is that C > 0. Assuming that < C < 1/2, 
we can determine the asymptotic behavior or infrared limit of the above J' functions, by 
taking the x — limit {p <^ A). We write 

Ji(C,x) MC)x-'-'^-B,iC) + C,iC)x'-''^ + ..., ^ = 1, 2, 3, (23) 

or 

JiC, x) ^ A{C)x-'-'^ - B{C) + C{Ox'-^^ + . . . , (24) 

with ^(C) = Eii A(C), B = Eii B,,C = ELi C„ and where 

11 1 
'^^ = 2(2TC)^ 2(2 + 0(1 + 20' -^'^ (2 + 0(1 + 20' ^ °' ^^^^ 
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and 

' F(-2-2C.-C;-i-2C;-i) rp + Ord + C) 



(2 + 20 ' ^' ^' ^' ' r(3 + 2C) 

r(-2-2C)r(i + c) l 
r(-i-C) J' 



(26) 



B.iC) - -^{(rio^(i + C,-C;2 + C;-i)-(^F(-i-2C,-C;-2C;^ 

r(i + c)r(i + c) r(-i-2C)r(i + c) l , . 

r(2 + 2C) ^ r(-c) J' ^ ^ 



and 



^3(0 = ^^^^|iF(C,-C;i + C;-i)-^i^(-2C,-C;i-2C;-i; 

_ r(C)r(i + c) r(-2C)r(i + c) l 
r(i + 2C) r(i-c) /■ 



(28) 



For < C < 1/2, only the leading scaling terms or power-laws and of Eq. ( p4|) are 

relevant for our purposes, therefore the functions Cj are not displayed. Then the vacuum 
polarization is 



^ / ^ 2A 



7r2 



1 / n \ 1+2C' 



1 + 2C ' ' VA 



(29) 



where 8 = 81 + 82 + 83. The critical coupling gc is given by the condition A^^(O) = in 
Eq. (0), thus 

- =Ai+A2 + A3 = (30) 

With the expression (|29|), the boson propagator (Eq. (^) takes the pure power-law form 
(0) at the critical coupling g = = I + 2(: 

consequently, the anomalous dimension rj = 1 — 2(. 

It can be shown that 83{( ^ 0) ~ and that 81 and 82 approach zero when C ^ 0, 

thus 8{0) = The case C = {i.e. A{p) = 1) in the above equation simply corresponds 

to the one-loop or leading correction, Ila{p) — (2A/7r^)(l — 8{0)p/A), see Ref. 0. 

At C = 1/2 it can be shown that 8 diverges, but so that 

8{C ^ 1/2) = C(C ^ 1/2), 8{C ^ 1/2) ^ (32) 



and J7'(C, x) reduces to 



:^(C,a;)^^ + ^lnx (33) 
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at ( = 1/2. Consequently, at ( = 1/2, the vacuum polarization Ho- of Eq. (H) reads 



^ / N 2A 



3 p' 



- + -^ln^ 



8A2 



V 



A 



(34) 



This means there is no pure power-law form at C = 1/2 and consequently the model becomes 
non-renormalizable in this case. This situation resembles the situation in 3 + 1 dimensional 
four-fermion models such as the gauged Nambu-Jona-Lasinio model, when the gauge cou- 
pling approaches zero, e.g., see Ref. ||18| . 

With the asymptotic form for the boson propagator, Eq. (|3ll) and the Ansatz (|TT1) for 
the wave function A, we can compute the Euclidean self-energy contribution Eq. (|10|) and 
verify the self-consistency of Eq. (pH]). Thus 

A 



A« 



4S(C) 



dk 



k ■ p 



47r p'^k-'^{k^ +p^-2k- p) V2+C • 



(35) 



Again the angular integral can be performed straightforwardly; 

A 



AC 



dk' 



k<- 







4S(C)2(3-2C)(1-2C) 

- + {i-2C)kp+p^]{\k~p\y-^'^ 

The integral is expressed as follows: 

^(C,P/A) 



k^ - (1 - 2C)kp + p^] [k + p)^-2C 



p 

8^(0(3-20(1 -20 VA 



(36) 



(37) 



where 



^(C,p/A) = J 



dk k^- 



p p 



p2 



;i-20- + i 
p 



1-2C 



1 + 



- + 1-20- + 1 



1-2C- 



dk k-<+'^ 



k p' 



-C+3 



1-2C 



P 



1-2C' 



(38) 



The function JF is written as 

1 1 
= J dtt^~^K{t) + j dtt^-^K{t), 

x 

where 



(39) 



K{t) = -{l-2C)t + l (1 + t) 



1-2C 



t^ + {l-2Qt+l 



t) 



1-2C 



(40) 
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Performing the integral, and determining the leading small x = p/A behavior of the function 
JF, we find for 

^ r(c)r(2-2C)(i + C)(2C-3) i ^^ 
^ - r(2-c) ^^0 — + ^^(C'-i + 2C,i + C,-i) 

(I - 2C) 1 

^ ^^-F(1 + C,-1 + 2C;2 + C;-1) + 7^F(2 + C,-1 + 2C;3 + C;-1) 



r(-i-c) C i-C 

fl — 2(^1 1 

+ F (-2 + C, -1 + 2C; -1 + C; -1) - j^F (-3 + C, -1 + 2C; -2 + C; -i) 

^^ 2(l-20(2C-^--5C^.2) ^^(^,,.) (41, 

Thus for small x = p/A, the can be written as 

^(C,x) ^ ^(C) + 0(x^), (42) 

where we write J-'iC) = J^{(,0). For values of ( with < ( < 1/2, the self energy has the 
scaling form 



where 



^ -A(C) (^^ ) , (43) 



8^(0(3-20(1 -20' ^^^^ 

Using the scaling form Eq. (^) and the Ansatz (|TTp in Eq. (^), we get the self-consistency 
relation 

f = 1 - NEa{p) = NXiO (1)"^ + 0{1). (45) 

This gives the eigenvalue equation (for < C < 1/2) 

A(0 = l/N. (46) 

This equation fixes the dependence of the anomalous dimension ( on A^. 
In the limit C ^ 0, the function JF given by Eq. (^) reduces to 

J^(0,x) ^ -21na;. (47) 

At C = 0, using B{0) = 7i^/8 and Eq. (g^), we find that S^(p) ^ 2/(37r2) ln(p/A) which is 
consistent with the one- loop or leading 1/A^ correction to S^i in the Gross-Neveu model [§. 
At = 1/2 the function JF(0 has a first order root and we can extract the residue by taking 
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0.2 0.4 0.6 

1/N 



0.8 



FIG. 3: Exact solution of Eq. ( ^q ) compared with = 1/2 — 5/(47rA^); the exponent C is plotted 
versus 



the limit C ^ 1/2 in Eq. (||) (at a; = 0) 



C-i/2 1 - 2C 



J dt [rl/2(l + i2)ij^i±^_ 2^1/2 





12-K 



t 

1+t 



+ Jdt r^/2(i + t2)ij,i±l_2t-5/2 



(48) 



Thus F(C) ^ 127r(l-2C)/5 for C T 1/2. Near C = 1/2, we obtain by using Eq. (ggD, Eq. 
and Eq. (H) that 



A(C) 



27r(l - 2C) _ 1 
5 ~ N 



2 AnN 



(49) 



The exact solution Eq. (|4q) and the above approximate solution are depicted in Fig. y. 
Clearly Fig. ^ shows that if N ^ oo then ^ — >■ 1/2. This behavior is a result of the fact that 
the functions 1/B{() as defined in Sec. [TV] and the function J-'{C) defined in Eq. (^) have a 
first order root at ( = 1/2. There, the eigenvalue equation is of the form given by Eq. (^^. 

Finally we mention that the leading scaling behavior is independent of the choice of 
momenta flowing in the self energy diagram Ha given by Eq. (0); an alternative choice. 



e.g., k + p for the fermion propagator and k for leads to the same scaling form (53) and 
eigenvalue equation (HBI). 
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V. THE GAP EQUATION AND THE CRITICAL EXPONENTS 



In the previous section the scahng laws for A{p) and Ao-(p) at the critical coupling g = Qc 
were obtained. This yielded the anomalous dimension rj = 1 — 2( with ( given by Eq. (^6[). 
In order to obtain the other critical exponents the set of SD equations should be studied 
either in the subcritical [g t gc) or supercritical {g J, gc) regime. In these regimes the 
analysis of the SD equations depicted in Fig. |1] is more complicated, due to the fact that a 
nonzero infrared mass scale (the inverse correlation length) enters the problem. For instance 
in the subcritical regime the inverse correlation length or is given by the complex pole 
p"^ ~ exp{—i6) in Ao-(p) |]20|. The critical exponent u can then be derived from the scaling 
law rricr ~ {gc — gY ■ If we consider Eqs. (|^) and (^) in the subcritical regime, it follows 
that V = 1/(1 + 2(). Consequently the anomalous dimension 7 = 1 as can be determined 
from the limit A^{p — ^0)~l/|5'c — with x the usual susceptibility P, [T^ 



In the subcritical regime the a boson propagator will be of the renormalizable form: 

^2 /A\^^^'' 1 



2Ai3(C) \pj [l + (m2/p2^(^+2^)/^ 

where rricr is the generic notation for inverse correlation length (rrio- oc rridyn ~ 1/^)- If this 
expression (^OD is inserted in Eq. ([I0|) a careful analysis will show that the fermion wave 
function A behaves as follows in the subcritical regime, 

(51) 

where c_ is a nonuniversal constant. This constant c_ should be determined from the 
nonlinear SD equations (0)-(0), but for the scaling arguments its precise value is irrelevant. 
In the supercritical regime the wave function will be of the same form as Eq. (^T]), but with 
a different constant of proportionality, e.g., c+ instead of c_. 

In rest of this section the supercritical regime [g > gc) is considered and we study the 
dynamical symmetry breaking SU(A^) — *-U(l). As was explained in Sec. |T|, in the leading 
large N approximation only the tad-pole graph is included in the gap equation for M. This 
gap equation with bare mass h reads (in Euclidean formulation) 

M = /i + 4G / — — . (52) 

Je (27r)^ A2(p)p2 + M2 ^ ^ 

For the wave function A the solution which is valid in the broken phase should be substituted. 
However, since the scaling laws for the order parameter and dynamical mass are mainly 
determined by the ultraviolet behavior mo- ^ p < A of the theory, we can safely use 
Eq. (0). Since the mass M in Eq. (|52D is independent of momentum, the above integral 
can be performed. Using Eq. (0) for the fermion wave function A, we obtain 

. h / A2 \ / 3 ^ 3 A2 \ 

The argument of the hypergeometric function approaches —00 as the critical curve is ap- 
proached {g ^ gc, M ^ 0) at h = 0. Therefore, close to criticality, the gap equation can be 
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approximated as follows: 



" '*T(i + .)r(i-.)f^r'-\ (54) 



\ Qc ) M 3' ' \ AJ 

where 6 = 3/(2 - 2C) and gc = l + K. 

Note that within our approximation M = (as). Therefore, treating M as the order 
parameter, the gap equation (|5^) is the quantum mechanical equivalent of the equation of 
state known in statistical mechanics and in analogy various critical exponents can be derived 
from it (see Refs. p|, |19|). For instance, the critical exponent 6 follows from the scaling law 
h^M^atg = g^; 

hocM^'-^+' =^ S = 2b-1='^ (55) 
and the critical exponent /5 follows from M ~ {g ~ gcY h = 0; 

g-g,r^M^'-^ =^ P = -^ = l^. (56) 
^ ^ 26-2 1 + 2C ^ ^ 

In a similar manner the critical exponent 7 can be computed from x = dM/dh\h=Q ~ 
{9 - 9c)~^; 

^ ^ M^-^' - g^)-^ =^ ^ = 1. (57) 

on h=o 

The other critical exponents can be derived from the effective potential and the scaling form 
for the composite boson propagator Eq. (31). Summarizing, the critical coupling gc and the 
critical exponents are expressed in terms of the exponent ( of the fermion wave function A 
which is determined by Eq. (^); 

= 1 + 2C, /? = (1 - 0/(1 + 20, S = {2 + 0/(1 - 0, 
z/=l/(l + 20, 7=1, ^ = 1-2C. (58) 

These exponents satisfy the three dimensional hyperscaling equations 7 = (3{6 — 1), 7 = 
z/(2 — rj), and Su = 2/3 + 7. The critical exponent u follows from the scaling law for the 
dynamical mass rridyn (inverse correlation length) describing the physical fermion mass in 
the broken phase close to criticality; m^yn ~ M^/^^~'^^ (see Sec. 0). 

In the limit N —* 00, keeping only terms to order 1/A^ by using the expression ( ^9]) for 
C, the critical exponents are 

These exponents clearly differ from the next-to-leading large exponents of the Gross-Neveu 
model 0; 



Note, for instance, the fractions of l/{7rN) appearing in the exponents (|59|), whereas the 
Gross-Neveu exponents contain fractions of l/(7r^iV). 
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FIG. 4: Next-to-leading order corrections to the Yukawa vertex 

From Eq. ( pT]) and (|3T|) , it follows that the anomalous dimension rj^ of the fermion field 
is 7]^ = — C/2 ~ — 1/4 + 5/(877 A^) and that the anomalous dimension rj^ of the composite 
boson field is r/o- = C ~ 1/2 — 5/ {Arc N). Since the /3-function of g vanishes at critical coupling 
9 = 9c {i-e., P{9c) = 0), the Yukawa vertex should satisfy the following Callan-Symanzik 
equation at g = gc- 

r]^ + 2r]^ + A j T^{k + p,k)^ 0. (61) 

Clearly, the bare Yukawa vertex approximation T^{p + q,p) = is consistent with this 
scaling constraint, since the scaling dimension of the bare vertex is zero in agreement with 
the equality rj^ = —2rj^. 



VI. DISCUSSION 

In Eq. (^), we have summarized the critical exponents obtained within our bare vertex 
approximation to order 1/A^. Together with the anomalous dimensions C, = 1/2 — 5/(47rA^) 
and = 1 — for respectively the fermion and a boson propagators for the model presented, 
this is the main result of this paper. However, the truncation scheme depicted in Fig. |I| 
describes merely the leading behavior in the 1/A^ expansion. Interestingly, the solutions of 
the truncated system depicted in Fig. |I] satisfy the hyperscaling equations for arbitrary but 
finite N . However, we cannot say with certainty that the values obtained in Eq. (^9]) are 
the correct critical exponents to order 1/A^, without a thorough investigation of the 0{1/N) 
diagrams for the Yukawa vertex. Examples of such 1/A^ Yukawa vertex corrections are 
depicted in Fig. ^. The factors of depicted in this figure represent index-loops which can 
be straightforwardly determined, using 't Hoofts double line representation (as explained in 
the Appendix). 



Suppressing factors of 1/A^ appear in the following way. From section using Eqs. (|3TD 
and (p^), we find that the a boson propagator at the critical coupling has the approximate 
form 

The eigenvalue equation (|49D identifies a factor of 1/A^ with 1 — 2^ (1 — 2^ oc 1/A^). Hence, 
each internal full a boson line in a Feynman diagram carries a suppressing factor 1/A^. This 
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illustrates that the diagrams in Fig. ^ are indeed of order Now the factor 1 — 2^ 

or is related to the Yukawa coupling gy mentioned in the Appendix. Since a four- 
fermion scattering amplitude, exchanging a virtual a boson would correspond to a factor 
gy, we can associate the square of the Yukawa coupling with gy ~ 1/iV, so that the 
factors gyN are indeed of order one, in retrospect. Using the above observation that each 
internal cr boson propagator corresponds to a factor and each index- loop to a factor A^, 
it is straightforward to verify that the vertex corrections depicted in Fig. ^ are order 1/N. 
Moreover, diagrams containing fermion loops with an odd number of out-going a lines are 
zero in the subcritical regime, due to restrictions imposed by the symmetry of the model. 

If = oo, we should get the exponents within the accuracy of the approximation. It 
turned out, however, that at A" = oo, ( = 1/2 and no pure power-law expression exists 
for the a boson propagator (see Eq. (p4|)). This logarithmic deviation from the power-law 
behavior causes the violation of hyperscaling and consequently spoils the renormalizability 
of the model. Such a violation could be a sign that the 1/A^ expansion is at the most 
an asymptotic expansion, since continuation to negative values of 1/A^ < are physically 
meaningless. But at the moment, this is a mere hypothesis. 

Nevertheless, when we take finite but large A^ then the model does have power-law renor- 
malizability and respects hyperscaling. The main question remains whether the exponents 
in Eq. ( pQ]) are correct to order 1/A^. In order to settle this question adequately one needs 
to investigate the 1/A^ corrections to the Yukawa vertex. We can think of three scenarios. 1. 
The next-to-leading 1/A^ corrections as depicted in Fig. ^ are finite and will contribute only 
to the "conformal invariant" or scale invariant part of the Yukawa vertex and do not alter 
its scaling dimension. This so-called scale invariant part of the Yukawa vertex depends only 
on the ratios of the in- and out-going momenta. 2. Each of the diagrams depicted in Fig. ^ 
is logarithmically divergent and the infinite sum of such diagrams give rise to a new scaling 
dimension of order 1/A^ for the Yukawa vertex. 3. The 1/A^ expansion breaks down due to 
the appearance of kinematic factors such as 1/(1 — 2^) terms, which give rise to additional 
factors of A^ and render the whole idea of topological 1/A^ expansion meaningless. In case 
that scenario 1 would turn out to be correct, then the set of critical exponents ( ^9l) will be 
correct to order 1/A^. Scenario 2 leads to an anomalous dimension for the Yukawa vertex of 
order 1/A^ and the expressions for the critical exponents and anomalous dimension are only 
valid in the leading order and up to logarithmic accuracy. In case that scenario 3 unfolds, 
our analytical approach has failed, since then we lost track of the perturbative parameter 
1/A^ and no standard techniques, except perhaps lattice simulations, are adequate to tackle 
the critical behavior of the presented model. By the way, each of three scenarios can be 
reconciled with the full conformal symmetry at the critical point of four-fermion models as 
advocated in Ref. 0. 

VII. SUMMARY AND CONCLUSION 

In this paper we have studied various aspects of dynamical symmetry breaking SU(A^) — >■ 
U(l) in a specific SU(A^) four-fermion model in the large A^ limit in 2+1 dimensions, by 
making use of 't Hoofts topological 1/A^ expansion. We advocated the apphcability of 
't Hoofts topological 1/A^ expansion to the presented model. Within the adopted 1/A^ 
expansion, the Yukawa vertex was approximated by its bare form. The main result of this 
paper is the determination of the exponent ( or anomalous dimension of the fermion wave 
function A{p) as a function of A^ in the planar approximation at the critical point g = gc- 
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The functional relationship between ( and given in Eq. (^6]) and Fig. ^. The approximate 
relationship was found to be C — 1/2 — 5/(47rA^). Subsequently the critical coupling gc 
and various critical exponents are given in terms of the exponent (. In the large limit, 
the critical exponents for this SU(iV) invariant four-fermion model differ considerably from 
the critical exponents of the Gross-Neveu model. The presence of A^^ — 1 light composite 
bosons turned out to drastically changes the scaling behavior of the fermion propagator in 
the vicinity of the critical point, giving rise to a new three dimensional universality class, 
see Eq. (^). Similar to the 2+1 dimensional Gross-Neveu model, the presented model is 
nonperturbatively renormalizable for large but finite N. The obtained anomalous dimensions 
of the fermion and a boson propagators are self-consistent with the bare Yukawa vertex 
approximation. 

For the future, it could be interesting to extend the presented approach to arbitrary 
space-time dimension d in order to investigate the perseverance of triviality at the upper 
critical dimension d = 3 + 1. 



Acknowledgments 



The author wishes to thank V.P. Gusynin for useful discussions and suggestions. 



APPENDIX: THE TOPOLOGICAL 1/N EXPANSION 



In this appendix, we illustrate the applicability of 't Hoofts topological 1/A^ expansion to 
our model Lagrangian (|l|). As is explicitly mentioned in Ref. ||ll|, the topological expansion 
is also valid for models, with a global U(A^) symmetry, having scalar fields with two U(A^) 
indices {i.e., the adjoint representation). Therefore, for large A^, we generalize the SU(A^) 
symmetry to U(A^). Subsequently, we perform a Hubbard-Stratonovich transformation, after 
which the Lagrangian (0) is expressed in terms of Hermitian scalar fields cTq^ = cr^^, 

C = tpaidtpa - i'a^pCTafB " -^CTapa/Sa, (A.l) 

with the flavor or "color" indices a, (3 running from 1 to A^. Contrary to the QCD and 
U(A^) gauge theories, "physical" observables are not restricted to be "color" blind. In other 
words, external fields or sources are allowed to carry U(A^) flavor indices. The Hermitian 
scalar field are linearly related to the A^^ real scalar fields (for U(A^)) of Eq. (j^). 



A=0 



with the U(l) generator = 5ap/\/N. 

Now the scalar fields a^/s are analogous to the gluon gauge fields of Ref. |ll|] , carrying two 
U(A^) indices, and consequently the topological 1/A^ expansion can be applied. Moreover, 
due to the absence cubic and quartic self interactions for the auxiliary scalar fields, the 
"flavor" Feynman rules are considerably simpler than the corresponding rules for the U(A^) 
gauge field theories. These flavor Feynman rules for the U(A^) symmetric phase are depicted 
in Fig. 1^, using the double line representation of 't Hooft, see also Ref. |T^. Each index 
loop in a Feynman diagram gives rise to factor A^. In Fig. |^, the propagators are fully 
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FIG. 5: The double index Feynman rules for the symmetric phase, with A^^{p) = —1/G + ^a{p) 
and S~^{p) = p — ^{p), where S is fermion self-energy. 

dressed. This is necessarily so, because of the fact that in a critical theory, with relevant 
four- fermion interactions close to or at a critical point, the vacuum polarization Il„ is of a 
comparable magnitude as the bare mass term 1/G. The fermion self-energy S is of the order 
of the canonical or free form p or larger (this will be shown in Sec. |V[). At the critical point 
the bare mass term l/G of the a boson propagator Ao-(p) is canceled by the momentum 
independent a boson vacuum polarization (11^) corrections and the remaining momentum 
dependent terms give rise to the anomalous scaling law for the scalar propagator. For the 



fermion propagator a similar argument holds, see also Sec. |IV. 

For d = A the leading expansion for U(A^) gauge theories consists of complicated 
planar networks of gluon exchanges. This is mainly due to the local nature of the symmetry, 
which introduces cubic and quartic gluon self-couplings and gauge-fixing interactions, whose 
particular structure depends on the space-time dimension. For d = 2, 't Hooft showed 
that the set of planar diagrams can be reduced to self-energy and ladder diagrams. As a 
consequence, the SDEs form a closed set and the d = 2 gauge theory is solvable in the 
leading expansion. 

Contrary to the local U(A'^) symmetries, the topological 1/N expansion is not space-time 
dependent for the global U(A^) symmetries of our interest, merely because of the absence 
of gauge fixing constraints or Faddeev-Popov ghosts. The space-time dimension d does not 
play a role in the selection rules for the planarity and thus the large N behavior of Feynman 
diagrams for global U(A^). For internal a boson exchanges is does not matter whether we use 
the real representation or the Hermitian representation for the composite fields to leading 
order in since 

A B A 

= iA„{p)6asS-yi3 + ^ i^lip) - i^aip)] Sap6^s, (A. 3) 

where the A^^ — 1 propagators are degenerate, = A^-, and A° is the propagator of 
the U(l) field ct°. In general A° 7^ A^., thus from Eq. ( |A.3|) it follows that the double line 
representation of the a boson propagator is violated at the most by a 1/A^ suppressed term. 

The planar expansion for the present 2+1 dimensional model is equivalent to the planar 
expansion for the d = 2 U(A^) gauge theory or 't Hooft model, since both models have a 
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iZ = 




FIG. 6: The planar approximation for the Yukawa vertex Tu, the a boson vacuum polarization Ho 
and the fermion self energy S. 





handle 



FIG. 7: The planar diagram corresponds to a factor A^(gyA^)'^. The hole diagram corresponds to 
a factor (gyN)'^, and the handle diagram to a factor l/N{g^N)^. 



single interaction vertex, see Fig. ^. This planarity can be verified order by order in the 
diagrammatic expansion of full a propagators A^-. We assume that with each internal full 
a boson exchange in a Feynman diagram an "effective" Yukawa coupling, gy is associated. 
We assume that combination QyN is fixed and of the order of one [|ry] , so that all diagrams 
with powers of QyN are of the same order and need to be summed nonperturbatively. The 
assumption that gyN is fixed for — * oo can only be validated in retrospect, since gy is not a 
free parameter like the gauge coupling g in QCD; gy is self-consistently determined from the 
planar SDEs and therefore depends on and G. For the a boson vacuum polarization Ho-, 
each planar diagram with n internal exchanges of a bosons is associated with a factor N"'g^. 
For the fermion self-energy E, each planar diagram with n a exchanges also corresponds to a 
factor N"'g^. Examples of vacuum diagrams with different topology (planar, hole, handle) 
are depicted in Fig. |^. 

For the fermion self energy, the planar expansion is generated by the rainbow approxima- 
tion. Moreover for the present model, a certain class of planar diagrams can be neglected. 
This class comprises all planar diagrams, which can be identified as vertex corrections. This 
can be understood in the following way. As depicted in Fig. ^ a full Yukawa vertex can either 
connect the left upper in-going line in 11^ with a right upper outgoing line, or connect the 
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upper left line with the lower left line. Whereas the bare vertex corresponds to the former, 
the ladder vertex corrections correspond to the latter. Since the generator r° is the only 
generator having a nonzero trace, it is straightforward to check that diagrams connecting 
upper left lines with lower left lines give only nonzero contributions to the propagator of 
the U(l) field, thereby causing the inequality A° ^ A„. For instance, in the real represen- 
tation, the vacuum polarization 11^ of the propagator gets a two-loop vertex correction 
which is proportional to J2b Tr This trace vanishes, except for A = 0. Using 

Eq. ( [A.3| ) in the fermion self energy, it follows that the contribution of A° is suppressed by 
a factor with respect to the degenerate SU(A^) propagator A^.. Therefore, we assume 

that all relevant planar diagrams, i.e., diagrams contributing to A^., are given by the bare 
vertex approximation depicted in Fig. ^. A similar planar diagram approach for the SU(A^) 
extrapolation of the Hubbard model was adopted by Foerster [O . 
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